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Abstract. In this article an intertwining operator is constructed which transforms 
the harmonic oscillator to the Dirac operator (the first order derivative operator). We give 
also the explicit solutions to the heat and wave equation associated to Dirac operator. As 
an application the heat and the wave kernels of the harmonic oscillator are computed. 
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1 Introduction 

The harmonic oscillators are interesting and important in their own right and play a 
fundamental role in the modeling of the quantum fields and are related to many mathe- 
matical and physical problems ([1], [2], [6], [7], [8], [9]). The aim of this paper is to give 
an intertwining operator T which relate the harmonic oscillator: 

i° = £-A 2 «>o (1-1) 

to the Dirac operator (first order derivative operator): 

D = ^ < L2 > 

We give also the explicit solutions to the following heat and wave equations associated to 
Dirac operator: 



(HD) 



(WD) 



d t U(t, X) = $?U(t, X) (t, X)eRlxM 
U(0,X) = U Q (X) U eC^(lR) 

d 2 t V(t, X) = £v(t, X) (t, X)eR* + xM 

V(0,X)=0 d t V(0,X) = Vo(X)eC^(M) 



Note that the Cauchy problem (HD) for the heat equation associated to the Dirac opera- 
tor D is nothing but the Cauchy problem for the transport equation. The wave equation 
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(WD) associated to the Dirac operator is considered in [3] p. 64 where the unequeness of 
solution is obtained but there are no known explicit solutions until now. 
As an application of our intertwining operator, from the heat and wave kernels of the 
Dirac operator, we give the explicit solutions of the following heat and wave equations 
associated to the harmonic oscillator: 



(HL a ) 



(WL a ) 



f d t u(t,x) = (jj^ - a 2 x 2 ) u(t,x) (t,x) e R* + x R 
\ u(0,x)=u (x) u eC^(M) 

dfv{t, x) = (j^ - a 2 x 2 ) v(t, x) (t, x) e M* + x M 

v(0,x) = d t v(0,x) = v (x),v e S a (M) 



with 



S a (R) = U : T 



e 2 



Note that the heat kernel for the harmonic oscillator is known for long time [1]: 
K a (x,x',t) = 



27r ^1 sinh(2at) 



exp 



--{x 2 + x u )coth(2at) 



sh(2at) 



(1.3) 



but to our knowledge the method used here and the obtained formula are new. 
For the wave kernel of the harmonic oscillator w a (t, x, x'), the following integral represen- 
tation is given for the wave kernel at the origin in [2]p.358: 



Wi(x, 0, t) = 



L 



4tt J c W 2zsinh(z) 



dz 



(1.4) 



where C is a contour symmetric with respect to the x -axis going throught the origin 
obtained by a smooth deformation of the circle C(\c, \c) of center c/2 and radius c/2 and 
for t < \x\ the kernel vanishes. 

Now we recall some facts about the Fourier transform: 
for / G L X (]R) the Fourier transform of / and its inverse: 



l 



(F- L f){x) = 



V2^ J- 
1 



f(x)dx 



^ x f{0di 



and we have the following formulas: 



a a 



](*) 



- x2 / 4s : s > 



(1.5) 
(1.6) 

(1.7) 
(1.8) 
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Recall also the complementary error function [3] p. 272 

Erfc(z) = / e- t2 dt 



and in terms of the Tricomi confluent hypergeometric U (a, c, z): 

Erfc(z) = ize- z V(f,3/2,z 2 ) = l e "* 2 £7(1/2, 1/2, z 2 ) 
2 The intertwining operator 



(1.9) 



(1.10) 



Definition 2.1 An operator T is said to be an intertwining operator if it relates operators, 
L and D, by 

TL = DT 



Lemma 2.2 For a > and <j> e S a (M) we have 

" 4>{x) = 



e 2 jr- 



Proof: Set 

we get 



e 2 



take the Fourier transform 

a 2 



_ii_i£2lli . „ J, i '°gisi 



0(x) =e" i/>(ar) 

a 2 



(5) 



(x) 



2 2 
— ax 



e 2 r/;(ar) = 



J" 



e 2 



9x 2 



2 2 
ax 



e°2 ^>(x) 



(0 



a 2 + 2ax TT + a 
ax^ ox 
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ip(x) 



-i z - 2a£— - a 



set 



we obtain 



e 4a "T 2 ^ 



e 2 



dx 2 



2 2 
ax 



e" tp(x) 



(0 = -2aC 



using (2.6)) we have 

-2<^f = 



il , ioalli 

;4a + 2 J 



dx 2 



2 2 
ax 



e—T- x \e 



>(0K*) 



Proposition 2.3 If / e 5(J?) and set 



(T0)(C)| C | =cxp (_2aX) = pJa + ! °2 £l JT 



e 2 



(2.1) 



(2.2) 

(2.3) 
(2.4) 



(2.5) 
(2-6) 



(2.7) 



(0 (2-8) 



(£)|£|=exp(-2aX) (2-9) 
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Then the operator T is an intertwining operator for the harmonic oscillator L a given in 
(1.1) and the Dirac opetrator D in (1.2): 

TL a = DT (2.10) 

Proof The formula (2.10) is a consequence of the lemma 2.2 

3 Heat and wave kernels for the Dirac operator 

Proposition 3.1 The Cauchy problem (HD) for the heat equation associated to the Dirac 
operator has the unique solution guiven by 

U(t,X) = U (t + X) (3.1) 

Proof: The Cauchy problem for the heat equation associted to the Dirac operator D 
reduced to the Cauchy problem for the classical transport eqution and its unique solution 
is given by (3.1). 

Proposition 3.2 The function 

w{t - x - x ' ] = 7m=^r (w^i) u (4 jp^jPf) (3 - 2) 

satisfies the wave equation associated Dirac operator D. 
and we have 

Wait, X, X') = ^=Erfc I , 1 ) (3.2)' 

V ^ \y/4:\X-X'\J 

where Erfc(x) is the complementary error function in (1.9). 

Proof: Set z = - 4{ /_ xl) we get f - - i{X -x>) &> 
ai . t 2 e 2 2 e _§_ t 2 d 

dt 2 4:(X-X>) 2 dz 2 i(X-X>)dz> dX i(X-X') 2 dz ' 

^X) - ^ {t ,X) = Y ^ Y \^ 2+ d/2 - = (3.3) 

the wave equation associated to the first order derivative operator is equivalent to the 
equation of confluent hypergeometric type [5] p. 268 

z<ft'{z) + (\/2-z)<ft{z) = (3.4) 

and an appropriate solution [5] p. 270 is z 1 / 2 exp (z)U (1,3/2, —z). 

Theorem 3.3 The Cauchy problem (WD) for the wave equation associated to the 
dirac operator has the unique solution given by: 

V(t,X)= [ W {X,X')V {X')dX' (3.5) 

J\X-X'\< | 
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where W{t,X,X') is given by (3.2) and (3.2)'. 

Proof In view of the proposition 3.2 it remain to show the limit conditions in (WD). 
And for this set X' = X + |s we have 

V(t, X) = OP* j\ -P (-^) U (l, 3/2, ^L) ds (3.6) 

where C = 2~ 1 / 2 -^|=. By the formula giving the first derivative of the Lommel confluent 
hypergeometric function [4] p. 265 

-^U(a,c,z) = -aU(a + l,c+ l,z) (3.7) 
az 

we obtain the limit conditions using the following behavior of the degenerate confluent 
hypergeometric function U (a, c, z), [5] p. 288 — 289 For z — > 0: 

[/(a, c, z) = (T(c - l)/r(a))js 1_c + 0(1), 1< 3£c < 2 (3.8) 

[/(a, c, z) = (T(c - l)/r(a))z 1 - c + 0(\zf c ~ 2 ), 3?c > 2, c ^ 2 (3.9) 

4 Heat and wave kernels for the Harmonic oscillator 

The Cauchy problem (HL a ) for heat equation associated to the Harmonic oscillator 
for u e L 2 (M) has the unique solution u belonging to C° ([0, oo[, L 2 (M)) see [6] more 
precisely there exists a semigroup (St)t>o of L 2 (K)-contractions such that for all t > 0, 
u(t, .) = StUo, the explicit solution given by 

/+oo 
K a (t, x, x')u (x')dx' 
-oo 

where the kernel K a (t,x,x') is given by the Mehler's formula (1.3) 

Theorem 4.1 If T is the intertwining operator given by (2.9) and if L a and D are 
the harmonic oscillator and the Dirac operator given respectively by (1.1) and (1.2) then 
we have: 

e tLa uo=T- 1 {e tD (Tu Q )} (4.1) 

5^„ = r->{2^ ( r„ 0) } (4.2, 

where e tA and s "^p' are respectively the heat and the wave kernels for the operator A. 

Proof The formula (4.1) is a consequence of the formulas (2.10) and the Cauchy problems 
(HD) and (HL a ). 

The formula (4.2) is a consequence of the formulas (2.10) and the Cauchy problems (WL a ) 
and (WD). 
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Theorem 4.2 The Cauchy problem (HL a ) for the heat equation associted to the 
harmonic oscillator has the unique solution given by: 

/+oo 

H a (t,x,x')u (x')dx' (4.3) 
-oo 

where 

H a (t,x,x') = 

^ - ^ at r 1/2 -p{ (e X e e-^ )2 + > 2 - x ' 2) ) (4 - 4) 

Proof: Using (4.1) of theorem 4.1 we have 

u(t,x) = e-"^ 12 / 2 /- 1 ^' 1 -"''^ 4 "/ (e- ax2 ' 2 u^ (te- 2at )](x) (4.5) 
u(t,x) = -J= e - at e ax2 / 2 x 

F- 1 [e- {1 - e ~ iat K 2 ' 4a ] * J- 1 [J" (e- ax2 ' 2 m ) (&- 2at )]{x) (4.6) 
using (1.7) and (1.8) we can write 



U(t,x) = ?—=-[ - V " — e-T^SF] * e 2at Uo(a;e 2a t)e -(a/2)*V- (4J) 

V 71 " V 1 - e~ 4a * 

= 6 ' V° / ^(x'e^^e-^^^e^^dx' (4.8) 
set ir" = x'e 2at in (4.8) we get the formula (4.3). 

Remark 4.3: The formula in (4.3) for the heat kernel associated to the harmonic 
ocillator agree with that given by (1.3) and the proof is left to the reader. 

Theorem 4.4 The Cauchy problem for the wave equation associated to the harmonic 
oscillator (WH) has the unique solution given by: 

v{t,x) = -^e™ 2 l 2 ^f^ I' Erfcf ^ V- 



J\lnt/('\<at \ V 2 ! In^/C'l / 

e«' 2 / 4a ,2, 

-^rHe™ /2 v*)(Z')d£,'](x) (4.9) 

Proof : The formula (4.8) follows from the theorem (4.1) and the Fubini theorem 
using the formulas (3.2) and (3.2)' and the fact that 

\Erfc(z)\ < ^ 
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We finish this section by the folowing corollary: 

Corollary 4.5 The heat and wave kernels for the cursin operator 

M = — + x 2 — 

dx 2 dy 2 

are given respectively by: 

i r+ca 

H a (t,x,y,x',y',t) = — / e^-y>H a (t,x,x')da 

w a (t,x,y,x',y',t) = — J e^-y^w a (t,x,x')da 

5 Directions for further studies 

We suggest here a certain number of open related problems connected to this paper, 
we are intersted in the heat and wave equations for the harmonic oscillator with an inverse 
square potential 

t HL ay f 9 t u(t, x) = (jj^ - a 2 x 2 - u(t, x) (t, x) £ R x R 
\ u(0, x) = u (x) u e Cq°(H) 



(WL a y 



dfu(t, x) = - a 2 x 2 - u(t, x) (t, x) e R* + x R 

u(Q,x) = d t u(0,x) = u (x) e C§°(R) 



Another possible extension is to consider the heat and wave equation associated to the 
power of the harmonic oscillator 



(HL a )' 



d t u(t, x) = [j^j - a 2 x 2 j u(t, x) (t, x) £ R* + x R 
u(0,x) = u (x) u £C§°(R) 



, WL ayi | d 2 u(t,x) = (j^ -a 2 x 2 ^ u(t,x) (t,x) £ R* + x R 

\ u(0,x) = d t u(0,x) = uo(x)£C^(R) 

Finally, we suggest a problems in direction of the non linear heat and wave equations for 
the harmonic oscillator and to look for global solution and a possible blow up in finite 
times. 

References 

[1]— Bcrlinc, N.,Getzler, E., Vergne, M. Heat kernels and dirac operator Springer Verlag 
2004. 

[2]— Greiner P. C. , Daniel Holocman, and Yakar Kannai, Wave kernels related to the 
second order operator, Duke Math.J.vol. 114, 329 - 387(2002). 



7 



[3]- Guelfand, I.M., and Chilov, C. E. Les distributions Tome 3 Theorie des equations 
diffcrcnticllcs, Dunod 1968 

[4]-Lebedev, N., N. ; Special Functions and their applications ; Dover Publications INC 
New York 1972. 

[5]-Magnus, W. Oberhittenger, F. and Soni R. P. Formulas and Theorems for the special 
functions of Mathematical physics, Springer- Ver log New- York 1966. 

[6]- Pazy, A. Semigroups of linear operators and applications to partial differential equa- 
tions Applied mathematical sciences 44 Springer Verlag 1983. 

[7]— Sheng-Ya Feng, A note on heat kernels of generalized Hermite operators, Tawanese 
Journal of Math. Vol 15,N°.5, 2035 - 2041(2011). 

[8]- Read ,S. and Simon B., Methods of modern mathematical physics II, Fourier analysis 
Sclf-adjointness, Academic Press, New York - London 1997. 

[9]-Thangavclu Hermite and Laguerre Semigroups some recent developpent Technical Re- 
port AT°2006/7, March 26,2006. 



Universite Gaston Berger de Saint-Louis B.P: 234. Senegal. 

E-mail adress: ahmcddou2011@yahoo.fr 

Universite de Nouakchott 

Faculte des sciences et techniques 

B.P: 5026, Nouakchott-Mauritanie. 

E-mail adresse: khames@univ-nkc.mr 



8 



